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Deterministic enhancement of the supercon¬ 
ducting (SC) critical temperature Tc is a long¬ 
standing goal in material science. One strat¬ 
egy is engineering a material at the nanometer 
scale such that quantum confinement strengthens 
the electron pairing, thus increasing the super¬ 
conducting energy gap A (iMgL: as was observed 
for individual nanoparticles [7|. A true phase- 
coherent SC condensate, however, can exist only 
on larg er s cales and requires a finite phase stiff¬ 
ness J fli| . In t he c ase of coupled aluminium 
(Al) nanograins (8l 4l0l |. Tc can exceed that of bulk 
A1 by a factor of three, but despite several pro¬ 
posals the relevant mechanism at play is not yet 
understood. Here we use optical spectroscopy 
on granular Al to disentangle the evolution of 
the fundamental SC energy scales, A and J, as 
a function of grain coupling. Starting from well- 
coupled arrays, A grows with progressive grain 
decoupling, causing the increasing of Tc- As the 
grain-coupling is further suppressed, A saturates 
while Tc decreases, concomitantly with a sharp 
decline of J. This crossover to a phase-driven SC 
transition is accompanied by an optical gap per¬ 
sisting above Tc- These findings identify granular 
Al as an ideal playground to test the basic mech¬ 
anisms that enhance superconductivity by nano¬ 
inhomogeneity. 

Bulk samples of pure Al represent a prototypical BCS 
superconductor (SC) with relatively low Tco ~ 1.2 K. 
Several studies since the late 1960s [iSi have shown 
a quite different situation for granular Al, i.e. thin films 
composed of 2 nm grains separated by thin insulating bar¬ 
riers, where a superconducting condensate is established 
via Josephson-coupling across the grain array. The cou¬ 
pling between the grains can be controlled during film 
growth, leading to samples with strong coupling and low 
resistivity (LR) in electrical transport compared to high 
resistivity (HR) samples with weak intergrain coupling. 
In LR samples Tc can be enhanced up to several times 
Tco, whereas it is suppressed to zero in HR samples, shap¬ 
ing a superconducting dome in the phase diagram, see 
Fig.IKa). 


To understand the behavior of Tc it is crucial to ac¬ 
cess the underlying SC energy scales associated with the 
amplitude and phase of the complex order parameter 
Ip = Ae"^. Indeed, while the SC energy gap A measures 
the pairing strength between the electrons, the true su¬ 
perfluid behavior can only be established if the Cooper 
pairs acquire the same macroscopic SC phase (p. The en¬ 
ergy scale controlling the rigidity of the condensate with 
respect to a deformation of this collective phase-coherent 
state is the so-called superfluid stiffness J. In ordinary 
BCS superconductors J exceeds A by orders of magni¬ 
tudes, and the SC transition at Tc is amplitude-driven. 
However, in the unconventional situation where A ex¬ 
ceeds J the transition is expected to be phase-driven, due 
to the loss of phase coherence at a temperature scale of 
order of J. Consequently, even though several finite-size 
effects have been proposed to explain the enhancement of 
A in isolated nano-grains [H-I^ , in agreement with recent 
observations in Sn nanoparticles Q, the behavior of the 
global Tc in large arrays of coupled grains is more difficult 
to assess, as here the grain coupling affects A and J in a 
competing fashion. On the one hand, in the LR regime 
the local pairing enhancement should be strong enough to 
overcome the smoothening of the local densi^ of states 
in the array due to strong grain coupling [6|. On the 
other hand as the film resistance increases, charging ef¬ 
fects will ultimately overcome the Josephson coupling be¬ 
tween grains, suppressing the global phase coherence. A 
direct implication of the above scenario is that HR gran¬ 
ular Al should undergo a direct superconductor-insulator 
transition (SIT), analogous to the one observed in homo¬ 
geneously disordered films of conventional superconduc¬ 
tors, like e.g. NbN, TiN and InOx il|,ll3|, or suggested to 
occur in unconventional cuprate superconductors, where 
J is suppressed by the proximity to the Mott-insulating 
phase |13l. Il4|. This analogy is made more interesting by 
the recent observation that in homogeneously disordered 
films and cuprate superconductors the SC 

properties become spatially inhomogeneous near the SIT. 
In this respect, the emergent granularity, either intrinsic 
or extrinsic, would constitute a general mechanism for 
the formation of superconductivity from an almost insu¬ 
lating normal state, irrespective of whether the driving 
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microscopic mechanism is based on disorder or correla¬ 
tion. In order to outline these potential analogies, and 
to better understand the mechanisms leading to a super¬ 
conducting dome in granular Al, it is crucial to assess 
experimentally the evolution of the characteristic SC en¬ 
ergy scales as a function of grain coupling, which so far 
has been only partly explored [ffl, [2l,[23]. 
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FIG. 1. (Color online) Superconducting dome and dy¬ 
namical conductivity, (a) Critical temperature Tc as a 
function of the normal-state resistivity (measured at 5K) of 
granular Al films studied in this work. Yellow symbols refer to 
the samples displayed in panels below. Tc encloses a dome-like 
superconducting regime with low-, optimal- and high resistiv¬ 
ity regimes, (b-d) (Normalized) spectra of o'i(i/) and 02 {v) of 
samples located on the left (sample 1), the right (sample 3), 
and at the maximum (sample 2) of the SC dome. The solid 
lines are fits to the Mattis-Bardeen theory. Note that the 
fit on (Ti disregards the low-frequency range due to excessive 
conductivity beyond Mattis-Bardeen theory. 


Thin films of 40 nm thickness were deposited on 
10 X 10 X 2 mm^ MgO substrates via thermal evap¬ 
oration. The degree of grain coupling was tuned by 
controlling the O 2 partial pressure during deposition 
and quantified by the dc-resistivity pdc in the normal 


state at 5K. Using a Mach-Zehnder interferometer we 
measured the complex transmission t of radiation in 
the spectral range 3 — 18cm“^ transmitted through the 
bilayer system and through a bare reference substrate, in 
order to disentangle Al and MgO contributions. Optical 
^He-cryostats allowed cooling the sample to T > 1.65K 
[ 2 ^. At the same time we measured the dc-resistance 
Rdc{T) (four-point geometry) to obtain Tc and to 
exclude heating upon THz absorption. In total, we 
examined nine samples with different resistivity values 
covering both sides of the SC dome, see Fig.[T](a). We 
measured t of all samples in the normal state well above 
Tc and in the SC state at Tbase = 1-65 K, as well as 
the temperature dependence of t for representative HR 
and LR samples. The complex dynamical conductivity, 
(t(i/) = +ia 2 {v) is calculated from the optical data 

using the Fresnel functions [H, . 


The great advantage of THz spectroscopy over other 
techniques is that it allows to simultaneously extract the 
fundamental energy scales of interest, A and J, from 
the measured cri{h') and 172 • Indeed, while the pair¬ 
breaking energy scale fixes the threshold for optical ab¬ 
sorption in <Ti{iy) below Tc, the superfluid stiffness J is 
directly connected to the inductive response (72 (j^)- In the 
present samples cri{v) and cr 2 {iy) can be adequately de¬ 
scribed by means of the Mattis-Bardeen (MB) equations 
for dirty superconductors. The accuracy of the MB fit of 
cri(z/) and (J 2 {v) is shown in Fig.[TJb-d) for representative 
samples covering the low- (sample 1), high- (sample 3), 
and optimal-resisivity (sample 2) regimes of the phase 
diagram. Even though the fit captures well the increase 
of conductivity at i/ > 2A/(he) (where h is the Planck 
constant and c is the speed of light), it underestimates 
the measured (7i(z/) at low frequencies. Such an excess 
conductivity strongly resembles the one observed, e.g., in 
disordered NbN and InO films Q and in cuprate films 
[i|, and is attributed to SC collective modes , not 
included in the MB theory. In the case of granular Al, 
where the Josephson coupling between grains is expected 
to be spatially inhomogeneous, this excess conductivity 
may be attributed to SC phase fluctuations, made opti¬ 
cally active by disorder |5l-l7|. The excess conductivity 
is seen in several samples at iz < 2A/(hc), irrespective 
of the grain coupling. A more quantitative study goes 
beyond the scope of this work and will be considered 
elsewhere. The gap value A(T) extracted from the fit 
of (7i(i/) is extrapolated to T = 0 assuming the BCS 
temperature dependence for A(T)/A(0), where the ra¬ 
tio A{0)/ksTc (with ks the Boltzmann constant) is not 
constrained to the weak-coupling value. 


The superfluid stiffness J is determined from the in¬ 
ductive response, i.e. (J 2 (v), which is proportional to 
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FIG. 2. (Color online) Superconducting energy scales. 
Tc, A(0), J(0), Ja( 0) (expressed in units of temperature) and 
A{0)/ksTc as a function of normal-state resistivity (measured 
at 5 K) of granular A1 films. Tc (green diamonds) encloses 
a superconducting dome with a maximum around 700 fiQcm 
where Tc is enhanced by nearly a factor of three compared to 
the bulk value. A(0) (olive stars) follows the increase of Tc 
on the left side of the dome for LR samples while it saturates 
in the HR regime. This is reflected in the ratio A{Q)/ksTc 
which increases from the weak-coupling value 1.78 (dotted 
line) to 2.25 when crossing from the left to the right side of the 
dome. The calculation of superfiuid stiffness from 02 ( 1 ^) and 
from A(0), i.e. J(0) and Ja( 0), is subject to a uncertainty 
reflected by the shaded area. J(0) follows approximately a 
1/pdc-behavior, as expected from Eq. (IS6I) for a constant value 
of A(0), and becomes comparable to A(0) in the HR regime. 
Dashed lines are guides to the eye. 


ris/m* [l^l36|. More specifically, we define 


h^nga 

4m* 


= 0.62 X 



( 1 ) 


Here, a is a transverse length scale, expressed in A, A is 
the penetration depth in and ng/m* = In 

an isotropic three-dimensional (3D) system, the length 
scale a in equation o is the SC coherence length ^ 0 , 
which is the natural cut-ofF for phase fluctuations, while it 
crosses over to the film thickness in the two-dimensional 
(2D) limit. Previous measurements of the upper criti¬ 


cal field in similar granular A1 samples [la . [Sggave an 
estimate of ^0 — 10 nm, while the analysis [3a, |38| of 
the paraconductivity above Tc indicates a 2D character 
with an effective 2D thickness for SC fluctuations of the 
order of ~ 15 nm throughout the phase diagram. As all 
these length scales are of the same order as the film thick¬ 
ness, we compare the evolution of the superfluid stiffness 
([1]) in our samples by choosing a constant value a = 10 
nm for the sake of simplicity. Once ng(T) is determined, 
the zero-temperature extrapolation follows from the two- 
fluid formula [ 3 ^ . 

The above analysis was applied to all samples under 
study. Fig. [2] comprises the results for the SC properties 
Tc, A(0), and J(0) (all expressed in units of temperature) 
as well as the ratio A(0 )/ksTc and presents them as func¬ 
tions of the normal-state resistivity. With increasing re¬ 
sistivity, Tc is first elevated from 2.7 K to a maximum Tc 
of 3.15 K before it is suppressed to 2.2 K. This SC dome 
with a maximum at about 700 ^Dcm is in good agree¬ 
ment with previous works on granular A1 composed of 
2 nm grains lH, 13 , 4^ . The enhancement of Tc in the 
LR regime is accompanied by a concomitant increase of 
A(0), so that the ratio A{0)/ksTc remains around the 
weak-coupling BCS value 1.78 for all LR samples. This 
behavior is in sharp contrast to the usual suppression of 
both Tc and A(0) for intermediate (but not too strong) 
disorder in homogeneously disordered films of conven¬ 
tional superconductors liB^. Thus, our measurements 
provide the first experimental confirmation that in gran¬ 
ular A1 the Tc enhancement is due to an increase of the 
local pairing scale A in each grain by finite-size effects, 
as suggested by several theoretical works in the past [l|- 
[^. However, as the dc-resistance increases further, phase 
fluctuations become more prominent and the overall Tc 
of the nanoparticle array is suppressed, even though a 
large local pairing survives. This is demonstrated in the 
upper panel of Fig. 2, where we show that on the right 
side of the SC dome A[Q)/kBTc increases up to around 
2.25, i.e., Tc is reduced more strongly than A(0), which 
remains fairly constant with increasing resistance. 

The crossover to a SC transition driven by the loss of 
phase coherence in the HR samples is confirmed by the 
strong suppression of T(0), shown in Fig. [51 The approx¬ 
imate scaling of J(0) with can be understood within 
the MB framework, where J(0) can be estimated from 
the A(0) obtained from ai{i/) and the normal-state re¬ 
sistivity as 


Ja(0) = 


Rc ^A(O) 
Rsa 4 ^ 


( 2 ) 


where Rc = hje^ and i?sq = Pdc/o, with same scale a 
as used in equation CD. Using A(0) extracted from tri 
we calculate a stiffness Ja( 0) that is similar to J(0) 
extracted from a 2 {r') following equation (HD , see Fig. [2] 
As the grain are progressively decoupled the 1/pdc 
prefactor in equation (IS6I) varies strongly, dominating 
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FIG. 3. (Color online) Temperature evolution of spectral 
gap. (a-b) Temperature dependence of normalized o-i(i/) 
of a granular A1 sample in the high- and optimal resistivity 
regimes (sample 2 and 3 in Fig.[T](d). In case of the HR sam¬ 
ple, the suppression of (Jiii') below Tc = 2.55 K (dashed lines) 
persists up to T = 2.8 K (solid lines), whereas the spectral gap 
closes right at Tc in the LR regime, (c) Temperature depen¬ 
dence of the spectral gap for samples from the optimal (stars, 
sample 2) and high resistivity regimes (diamonds, sample 3). 
The blue data traces pdciT) of the HR sample. For the HR 
sample, deviations from the BCS prediction for A(T)/A(0) 
(black solid line) appear already at T/Tc < 1, where A is 
anomalously large. The persistence of a gap across Tc (empty 
diamonds) is in striking resemblance with strongly disordered 
or correlated superconductors. 


the overall scaling of J(0). As a consequence, even in 
the LR regime the absolute value of J(0) is considerably 
lower than what is expected for conventional clean 
superconductors, where J(0) scales with the Fermi 
energy. According to equation o this translates to 
a penetration depth for LR samples of A ~ 1 /rm 13611 
much larger than A ~ 50 nm reported for bulk A1 [4I| . 
However, while in LR samples J(0) still exceeds A(0) 
by one order of magnitude, in HR samples J(0) is 
strongly suppressed and becomes comparable to Tc- 


We notice that the given estimate of J(0) should be 
taken as an upper bound, since it neglects the additional 
reduction due to inhomogeneous phase fluctuations 
0. Howeve^ the comparison with previous SQUID 
measurements [ 3 , of the penetration depth suggests 
that this effect is still quantitatively small for the sam¬ 
ples under consideration. The slightly different values 
of J(0) and Ja( 0) is reflected in the shaded area in Fig. [21 


Apart from the anomalously large value of A(0)/fcBTc, 
in the HR regime we find further evidence for an un¬ 
conventional behavior of superconductivity from the 
dissipative conductivity. Fig. m a) and (b) compare 
(normalized) cri(z/) spectra of high- and optimal resis¬ 
tivity samples (sample 2 and 3 in Fig.[T]) at various 
temperatures below (dashed lines) and above (solid 
lines) Tc- In the HR regime, MB theor y ag rees with 
the measured data at low temperatures [36|, whereas 
by approaching the data deviate significantly from 
the MB prediction, with a strong suppression of cri(z 2 ) 
at low frequencies, see Fig.|3Ka). This suppression 
in <Ji{v) exists up to T = 2.8 K well above the SC 
transition, as evident from Fig. EDc) where Pdc{T) of 
the HR sample is shown for comparison. In contrast, 
in the case of the optimal-resistivity sample (Fig.|5jD) 
both <Ti{v) and i spectra [s^ contain no signs of a 
spectral gap above Tc, and A(T) follows closely the BCS 
temperature dependence up to Tc, see Fig.|3t. For the 
HR sample the quality of the MB fit degrades already 
at T/Tc — 0.8, as signalled by the larger error bars in 
A(T) reported in Fig.jS}:, and A(T) evolves smoothly in 
a finite spectral gap found up to the highest measured 
temperature. The same anomalies are observed in the 
analysis of the paraconductivity, discussed in details in 
the Supplementary section 36|. In particular, while for 


low- and optimal-resistivity sample the paraconductivity 
is well described by ordinary Aslamazov-Larkin type 
of Gaussian SC fluctuations, in the HR samples we 
observe clear deviations that are in agreement with 
previous observations of an unconventional fluctuation 
regime based on magnetotransport and Nernst effect 
[l3l [. All these findings suggest that the anomalous T 
dependence of A and cri(i/) for HR samples near Tc 
can be attributed to a pseudogap above Tc, which can 
be viewed as a direct and natural consequence of the 
phase-driven transition at high resistivities, even though 
a full theoretical understand of it is still lacking. 


From our comprehensive measurements of the dynam¬ 
ical conductivity of superconducting granular A1 thin 
films at THz frequencies we determined the dependence 
of the energy scales Tc, A(0), and J(0) on the decoupling 
of the A1 grains. We show that decoupling promotes 
the individual nature of the grains and enhances the 
local pairing amplitude in each grain due to finite-size 
effects. The enhancement of both Tc and A(0) in the 
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low-resistivity regime is eventually overcompensated in 
high-resistivity samples by enhanced phase fluctuations, 
which suppress while the pairing amplitude remains 
large. The strong suppression of J(0) and the persistence 
of a spectral gap above Tc in the high-resistivity regime 
indicate a crossover to a phase-driven transition. Our 
results provide clear constraint on further theoretical 
modelling of the superconductivty in granular arrays, 
that is a prerequisite to achieve full control in the 
Tc enhancement via engineered inhomogeneities in 
low-dimensional superconducting nanostructures. 
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Analysis of the conductivity and determination of A(0) and J(0) 


We analyze both and fT 2 (*^) by means of the ordinary Mattis-Bardeen (MB) formulas [ij for the optical 

conductivity: 




TTV 2 7 °° 

= ^5{v) + — dEg{E)[f{E)-f{E + hu)]- 
m*(7n hv 


Q{hv — 2A) 
hv 

E{E + hv) + 




hv 


[ dEg{E) [1 - 2J{E + hv)] (SI) 

•/ A — hiy 

(S2) 


where the function g{E) is 


9{E) 


E{E + hv) + A2 

VE^ - AV(^ + ' 


(S3) 


Here (7„ indicates the normal-state conductivity, that is featureless in the THz frequency range for our samples. 
To exclude any residual fluctuation effect we take as (t„ the one measured at 5 K, well above for all samples. 
In Fig. [ST] we show the accuracy of the fit of cri(j^) and a 2 {v) of a representative high-resistivity (HR) sample at 
several temperatures. We remark that the fits show deviations at frequency v < 2A/(/ic), where we observe an 
excess conductivity analogous to the one observed in other conventional superconductors as NbN and InO [3, |^. 
This has been attributed to superconducting collective modes not included in the MB theory, whose detailed 

study goes beyond the scope of this work and will be considered elsewhere. The obtained values of A(T) are fitted 
with a conventional BCS temperature dependence for A(T)/A(0) vs T/Tc, where IS.{{))/ksTc enters as the only free 
parameter. The BCS expression reproduces accurately the behavior of A(T) for samples on the left side of the 
superconducting dome, while for HR samples some deviations occur near Tc, that can be attributed pseudogap above 
Tc in HR samples. 

As a check of consistency of this procedure we remark that the emergence of a global superconducting state can 
directly be seen from the raw transmittance (transmission amplitude) |t| of the electromagnetic radiation passing 
through the A1 film. Below Tc and in the spectral range relevant for this work, \t\ is higher than the normal-state 
transmittance |to| and frequency dependent with a maximum close to the spectral gap 2A, i.e. the normalized 
transmittance t = |t|/|to| is above unity. The presence of a pseudogap above Tc manifests as t = |t|/|to| > 1 in the low 
frequency limit, or, equivalently t = 1 in case of no pseudogap present. In Fig. 3 we show t for a low-resistivity (LR) 
sample (pdc = 516/xHcm) on the left side of the superconducting dome for various temperatures above and below 
Tc = 3.15 K. While at low temperatures t is well above unity, any signature of superconducivity vanishes directly 
at Tc- This is contrasted by a HR sample. Fig. 3b, where we find t exceeding unity already at 2.8 K well above 
Tc = 2.55 K which suggests the opening of a pseudogap that smoothly develops into the superconducting gap. 

The superfluid density is given by the missing finite-frequency spectral weight between normal- and supercon¬ 
ducting state ai{v > 0), and it controls the weight of the zero-frequency delta response in Eq. (ISTI) . Since it cannot 
be directly accessed by optical probes, we extract it from the Kramers-Kronig transform fT 2 (*^) as 


Us 


27rm* 


lim v<72M, 


(S4) 


where m* and e are the effective electron mass and charge. We fit the experimental data in the THz range to the MB 
formula (IS2D . see Fig. ISH and extrapolate fit times frequency to zero. The value at zero temperature is then obtained 
from a standard two-fluids temperature dependence: 


ris 


1 - 



(S5) 


where rig^o is the zero-temperature superfluid density. As in the case of A(T), we find representative samples on 
both sides of the superconducting dome following this universal dependence on T/Tc, see Fig. IS3b b which justifies 
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FIG. SI. (Color online) cri(i/) (blue triangles) and a 2 {v) (green triangles) of a HR sample as functions of frequency at various 
temperatures in the superconducting state together with fits to MB formula (EU and (IS2II . Apart from moderate excess 
conductivity at subgap frequencies (denoted by the dashed line, excluded from the cri fit) MB theory reproduces our results 
remarkably well. Corresponding spectra of an LR sample can be found elsewhere @]. 


to calculate rig^o from Eq. (IS5I) for samples where no full T-dependence was measured. We remark that even though 
the exact value of m* is needed to extract Ug from the measured cr 2 {v), it is instead irrelevant when computing the 
superfluid stiffness J(0) = fi^ns^oa/Am*, which only depends on the ratio ns/m*. We also notice that the separate 
analysis of < 72 ( 1 ^) to obtain ng^ is in principle redundant, since within the MB theory we could use directly the value 
of A(0) extracted from the cti fits to determine the zero-temperature inductive response. Indeed, from Eq. (IS2I) one 
immediately sees that 172 ( 1 ^ 0,T = 0) = 7 rA( 0 )t 7 „/(/ii^). Since in the dirty limit t 7 „ coincides with 1/pdc in the THz 

frequency range, we can estimate Ug^ as: 


^ 27rm* 7rA(0) 


so that the corresponding estimate Ja of the stiffness corresponds to Eq. (2) of the main text. However, since the 
deviations of (7i(i^) from the MB behavior occurs exactly below 2A, we analyzed f72(i^) independently on f7i(zz), and 
we cross-check afterwards the consistency between the two approaches, see Fig. [S3k). Here the extracted values of 
ns,o and nf^Q are converted to the penetration depth A = poe^Us/m* (with jm the vacuum permeability), in order 
to compare them with direct measurements of A done in previous works [^, [l^ . Both estimates of A are consistent 
with each other and they are in very good agreement with previous findings. This shows also that the quantitative 
suppression of ng^ due to the collective-mode contribution below 2A in fJ 2 {v) is relatively small, and it justifies the 
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FIG. S2. (Color online) Normalized transmittance t = |t|/|to| of granular A1 films versus frequency and selected temperatures. 
The thick solid line sketches the temperature dependence of the dc-sheet resistivity pdc-{a.) Low-resistivity {Tc = 3.15K,pdc = 
516pflcm) and (b) high-resistivity sample (Tc = 2.55K, pdc = 2107/rr2cm). For the LR sample, we find |t|/|to| to increase 
above unity only below Tc in sharp contrast to the HR regime, where \t/to\ > 1 already above Tc indicates a pseudogap. 


use of the MB formula to extrapolate cr 2 to zero frequency. 


Analysis of the paraconductivity above Tc 


To analyze the effect of the superconducting fluctuations above Tc (here defined as the temperature where the 
resistance becomes immeasurably small) we focus on the Aslamazov-Larkin (AL) contribution to the paraconductivity 
Act [l^ . that is the most relevant near Tc- In agreement with previous w ork|l3 . we found that in our samples the 
paraconductivity has the temperature dependence expected in 2D, so that |12l| : 


Actal 

CTtv 
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_ 1 
IQhdALf^N e 16 I?c e ’ 


(S7) 


where Rc = = Pdc/dAL and dAL is a transverse length scale, of the order of the film thickness, that 

determines the effective 2D unit for SC fluctuations in our granular samples. The parameter e contains the temperature 
dependence, and in the BCS limit it is given by 



(S8) 
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FIG. S3. (Color online) (a) Penetration depth A versus normal-state resistance. The present work (colored dots) very well 
reproduces results obtained in previous works by of Abraham[^ and Gershenson [lo)| (empty and filled squares), where the 
inverse penetration depth has been directly measured. In addition, we note that the absolute numbers of A obtained from 
the inductive response (blue dots) are nearly identical with the calculation of A from A (green dots) within MB theory, (b) 
Superfluid density Ua of a LR and HR sample versus temperature. The solid lines are fits to the two-fluids approximation (IS5II 
to obtain the zero-temperature extrapolation riafi. For the LR and HR samples we estimated a mass enhancement of m*I'm = 2 
and 5, respectively, from the normal-state resistivity [ll| . 


We fitted all the data using the 2D AL formula (1571) . using dAL and as a free parameters. As one can see in Fig. 
[Sllthe agreement with the data in the LR regime is remarkably good, with an effective thickness dAL — 14 — 17 nm 
as an adjustable fit parameter, that is within a factor of 2-3 from the real film thickness. Given the granular nature 
of the film this is a reasonable approximation, considering that the fit reproduces the data up to temperatures as 
large as twice without any other adjustable parameter. On the other hand, when one analyzes the HR films, two 
remarkable differences arise: (i) the resistivity is not completely saturated up to temperatures as large as twice Tc, (ii) 
the fit with the AL formula fails around T ~ I.2Tc, since the experimentally measured paraconductivity decays faster 
than predicted by Eq. (IS7I1 . Interestingly, the very same behavior has been observed also in underdoped cuprates for 
samples in the pseudogap regime. In this case, it has been observed that a very good formula that works both near 
and far from Tc is the following one 


14|: 


^ctal _ Rg _ 1 ^ 

aN 16i?c eo sinh(e/eo) ’ ^ Tc’ 


(S9) 


that reduces to the usual one dHil when e <C eo, so that sinh(e/eo) — e/eo, but decays faster for e ^ eo- As it has 
been discussed in the context of cuprates [l3. [l5| . such a suppression of paraconductivity with respect to the standard 
formula (IS7I) can be indeed explained assuming that a pseudogap survives in the electronic Green’s function up to 
a temperature T* larger than Tc- In the case of granular A1 the formula (jS9() works remarkably well for the two 
most disordered hlms, see Fig. [SH where pseudogap signatures have been reported also from the MB analysis of the 
optical spectra. 
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FIG. S4. Fit of the resistivity with the AL formula (IS7I) for LR samples and with the modified AL formula (IS9II for the HR 
samples. The parameters of the fit are shown in each panel. 
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